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Abstract 



We solve the two-impurity two-channel Kondo model using a combination of 
conformal invariance and bosonisation techniques. The odd-even symmetric 
case is analysed in detail. The RKKY interaction turns out to be exactly 
marginal, resulting in a line of non-Fermi liquid fixed points. Explicit formu- 
lae are given for the critical exponents and for the finite-size spectrum, which 
depend continuously on a single parameter. The marginal line spans a range 
of values of the RKKY coupling / which goes from the infinitely strong fer- 
romagnetic point I = —00 (associated with a 4-channel spin-1 Kondo model) 
to a finite antiferromagnetic critical value I c > beyond which a Fermi liquid 
is recovered. We also find that, when the odd-even symmetry is broken, the 
marginal line is unstable for ferromagnetic /, while for antiferromagnetic / it 
extends into a manifold of fixed points. 
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The effect of inter-impurity interactions on quantum single-impurity models possessing 
a non-Fermi liquid ground-state is of crucial importance for the possible experimental real- 
izations of such systems [[J . It is also of considerable theoretical interest for understanding 
non-Fermi liquid behaviour in lattice models of correlated fermions starting from a local 
point of view |2||| . 

The two-impurity Kondo model with two channels of conduction electrons is one of the 
simplest model where this problem can be addressed. For a single impurity, this model is 
controlled by a non-trivial fixed point [|j], resulting in a specific heat coefficient C /T and 
susceptibility Xim P diverging logarithmically as T — > 0, and a universal finite-size spectrum of 
excitation energies differing from the free-fermion form ||. The corresponding two-impurity 
model is the simplest situation which brings in the competition between the formation of 
this non-trivial Kondo state and the ordering of the impurities via the RKKY interaction. 
It has been recently studied by numerical renormalization-group methods (NRG) f6^f§. 

In this letter, we present an analytic solution of the low-energy universal properties of 
this model using a combination of conformal field theory |5||| and bosonisation methods 
njnj- We find that the RKKY interaction (as well as other interimpurity couplings) is a 



marginal perturbation, giving rise to a continuous family of non-Fermi liquid fixed points. 
The finite-size spectrum and the critical properties vary continuously with the strength of 
the interaction. We obtain analytic formulae for this dependence. These results are in 
excellent agreement with recent NRG results They should be contrasted with the 

single-channel case in which, in presence of particle-hole symmetry, Kondo screening always 
dominate over RKKY ordering or vice- versa, resulting in two stable Fermi liquid fixed points 
separated by an unstable non-trivial critical point [0.0] • 

We formulate the model in terms of left-moving (chiral) fermions ipu a {x) on the full axis 
— oo < x < +00. I = 1, 2 is an index labelling the two impurity sites, % = 1, 2 is a channel 
index and a a spin index. The hamiltonian is written as: 

H = iv F Ezi Q dxipj ia (x)£ipu a (x) (1) 
+ J + (Si + S 2 ).(J 1 (0) + J 2 (0)) (2) 
+ J m 0i-S 2 ).(J 1 (O)-J 2 (O)) (3) 

+ + s 2 ). Zi,aMia(Q)¥^(ty + V>L(o)%MM°)) (4) 

+ (5) 

In these formulas, J\[x) = J2i,a/3 4 ! iia( x )~^f~' l f ) iii3( x ) denotes the spin-current at position x for 
species I. Our notations follow closely those of Ref. |J. Alternatively, one could work in the 
even/odd basis ip £j0 with respect to the midpoint between impurities. The combinations ipip 
correspond to (ip e ± ip a ) / v^2, respectively, and a parity transformation amounts to exchange 
the indices I = 1,2 for both impurity spins and conduction electrons. In order to make 
contact with the couplings Y eo>m used in Ref. 0, let us note the identifications: J m oc T m , 

j + oc r e + r Q , j_ oc r e - r Q . 

We shall start by identifying the global symmetries of the hamiltonian. For most of this 
paper we shall concentrate on the case J_ = 0, corresponding to a hamiltonian invariant 
under odd-even exchange (T e = T Q ). H has a higher symmetry in that case, with independent 
charge and channel (or 'flavour') transformations allowed for I = 1, 2: 
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ipiia -> e iei ijj Ua (6) 
hia^EjUg^a, U®eSU(2) (7) 

When J_ 7^ 0, only Q\ = 9 2 and = are allowed. The symmetry of the spin sector 
is conveniently discussed by looking first at the case of two decoupled impurities (each one 
interacting with two conduction channels) obtained by setting 1 = and J + = J m in addition 
to J_ = 0. In that case, H has independent spin-rotation symmetry: 

^ua - E kSVw , sf - £ w )sf (8) 

,a b 

where R ab (V) = l/2tr(a a Va b V^) , (a, b = x,y,z) is the adjoint representation of V G 
SU(2). Hence, two decoupled 2-channel Kondo models have global symmetry (SU(2) spin ® 
SU(2) flavour <£> U{l) c harge) 2 ■ Coupling the two impurities (I 7^ 0, J+ 7^ J m ) while keeping 
J_ = leaves unchanged the independent charge and flavour symmetries, but reduces the 
spin symmetry to the diagonal SU{2) corresponding to V^ 1 -* = in Eq.(^). 

At a fixed point, these global symmetries are promoted to local conformal symmetries 
13| . For decoupled impurities, the symmetry algebra consists imtwo copies of a product of 



Kac-Moody algebra for spin, channel and charge: (SU 2 (2) S ® SU 2 (2) f ® U{l) c ) 2 . (SU k {2) 
stands for the level-fc SU(2) Kac-Moody algebra, corresponding to the commutation relations 
of the sum of k independent SU(2) currents). When coupling the impurities with J_ = 0, the 
diagonal 577(2) symmetry of the spin sector gives rise to a S77 4 (2) algebra. The generators 
of this algebra are the sum of the generators of the two SJJ 2(2)^ for each impurity, that is 
the sum of the spin currents Ji(x) + ^(x). Hence, we must understand how the product 
SU 2 (2) S g) SU 2(2)3 can be decomposed into SU 4: (2) S plus some residual degrees of freedom. 
The answer is given by the so-called coset construction |L4] : 



SU 2 (2) S ®SU 2 (2) S = SU 4 (2) S ® A(2,2) (9) 
The algebra A(2, 2) turns out to be a N = 1 superconformal unitary model [T^llll corre- 



sponding to the m = 4 member of the discrete series with central charge c = |(1 — m ^+2) )' 
and thus has c = 1. This construction generalizes to the two-channel case the one made by 
Affleck and Ludwig in their solution of the one-channel two impurity problem [[|. There, 
the coset construction is S77i(2) s <g> SUi(2) s = SU 2 (2) S ® A(l, 1), where the algebra A(l, 1) 
is actually an Ising model with c = 1/2. 

This coset construction can be understood more explicitly when dealing with spin cur- 



rents (i.e for the adjoint representation of the algebra). Let us first recall [13| that the 
SU 2 (2) spin current Ji"(x) (a = x,y,z) for a given I = 1,2 can be represented in terms of 
three Majorana (i.e real) fermions x{ v ' z as follows: 

J?(x) = te abcX b Xi (10) 

This is particularly transparent when using the Emery-Kivelson bosonisation approach to 
the two-channel Kondo model 0,0, in which case \ x — sin$ s ,x y = cos§ s ,x z = cos$ s / 
where $ s ,$ s / are the boson field introduced in ref. [IIJ corresponding to spin and 'spin- 



flavour' degrees of freedom. Here, we are dealing with two copies of SU 2 (2) and hence with 
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six Majorana fermions. We combine them into three Dirac fermions and bosonize these new 
degrees of freedom as: 



X?(x)+^(x) = e^ ,a = x,y,z (11) 

Our notations are such that the free boson correlator reads: < $(r)$(0) >= —ln(r), so that 
e lfc * has dimension k 2 /2. In terms of these fields, the total spin current corresponding to the 
diagonal S77 4 (2) algebra reads: J x = + J% = cos(§ y — $> z ) (and cyclic permutations). 
It is convenient to introduce three linear combinations of boson fields as follows: 

$ = 73 ($x + ® y + $z) ' ^ = 72 ($x " ® y) ' v = 7e^ x + $y " 2$z) (12) 

In term of these combinations, the components of the total spin current read: 

J* = COS ^~^ U ) , jy = COs( fl+ ^ U ) , J' = C03(V2fi) (13) 

v 2 v 2 

Note that $ does not enter these expressions. Hence the two bosons /i, v are sufficient to 
describe the SU^iT) algebra (as expected from its central charge c = 2) and <3> corresponds 
to the residual v4(2,2) degree of freedom (c = 1). (The central charge c = 3/2 + 3/2 has 
thus been distributed as c = 2 + 1 in the coset construction Eq.((|). 

Thus, a very useful explicit realization of the algebra A(2, 2) as a free field theory of a 
single compact chiral boson $ has been found. Since $ x ,y,z have periodicity 2ir, the radius 
of $ is found from eq. (0) to be R = V3, which means that $ and $ + 2%R are identified. 
However, this construction based on currents in the adjoint representation does not reveal 
the full structure of the SU2(2) and A(2,2) algebra. Indeed, the three Majorana fermions 
involved in the currents of STJ 2(2) are not fully independent. For example, if we change their 
boundary conditions from periodic to antiperiodic, all three must be changed simultaneously 
so that the boundary condition on the current remain periodic. The boundary conditions 
on the two sets of Majorana fermions xf> I = 1,2 are independent however, which means 
that the fermions constructed in Eq.flTTD are not strictly speaking Dirac fermions but rather 
that $ and — $ must be identified. A(2, 2) is thus really a (chiral) 'orbifold' theory. This 
is of crucial importance in determining the operator content of this algebra, and thus the 
finite size spectrum of the present model. It implies that, in addition to the operators 
e ±«(nv / 3+m/2 v / 3)* anc } QnQ (^ w ith n,m integers), A(2, 2) contains two operators of dimension 
1/16 which do not have a simple boson representation. These operators are analogous to the 
(twist) operator associated with the order parameter for the case of the Ising model, which 
change the boundary condition of the Majorana fermion from periodic to antiperiodic flIS 



For convenience, the full set of primary operators of the A(2, 2) algebra and their boson 
representation (when it exists) is given in Table 1. 

We now consider the effect of turning on the 'RKKY' coupling /, and/or setting J + ^ 
J m , starting from decoupled impurities (but keeping J_ = 0, i.e T e = T ). An order by 
order perturbative calculation in these couplings (e.g of the free energy) involves correlation 
functions of the operator S±. S2 which is a product of spin correlators at the decoupled 
impurities fixed point. As will be detailed in a longer publication |I7| |, it can be shown 



that, for the purpose of calculating these correlations, the impurity spin Si (resp. S 2 ) 
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can be replaced by aiXi (resp. CI2X2) , where 01,2 are local real fermions needed to ensure 
proper commutations. Hence, the perturbing term of lower dimension associated with the 
RKKY interaction reads / dtaxa^Xi-Xz- m the bosonic language above, this translates into 
an induced boundary term in the A(2, 2) sector of the hamiltonian: 

Ham = V ~fJ dx(U(x) - + I{d+d - 1)^(0) (14) 

where we have set d + = (a t + ia 2 ) /y/2, n is the field conjugate to $ and I is some (non- 
universal) function of I and J + — J m . Hence the RKKY coupling is associated with a 
dimension 1 operator and is an exactly marginal perturbation. 

The hamiltonian ( |i~4"| ) is very similar to the X-ray edge hamiltonian (in presence of the 



charged core) in the bosonised form p8| . In the X-ray edge problem, the boundary term 
gives a phase-shift, which changes the dimensions of the operators e ±ifc * from y to ^(/c±^) 2 . 
This can be obtained by formally fusing e ±jfc * with Us = exp(i^-(d + d — l/2)$(0)). This is 
the analogue of the "fusion hypothesis" introduced by Affleck and Ludwig for the Kondo 
problem This fusion rule can be generalised to the present problem (with 5 some 

non-universal function of / and J + — J m ) provided we know how to generalise it for the 



operators of dimension (^). It turns out |L7j that the fusion gives back the (^) operators 



plus operators of dimension ^ + ~ and conformal descendants of these two operators. 

Using this method, we have derived the finite-size spectrum of the model for / / 0, 
J+ 7^ Jm (but J_ = 0) as a function of the single parameter 5, starting from the (known 
||) finite-size spectrum of two decoupled two channel Kondo models. This is displayed in 
Table 2 (where the parameter x = 2y / 3^ has been used). Note that the ground-state is the 
triplet of lowest energy for x > (ferromagnetic coupling), and the singlet of lowest energy 
for x < (antiferromagnetic coupling). Accordingly, the normalised excitation energy of a 
given state, (with L the radial length of the bulk system and Vp the Fermi velocity), is 
obtained from the total dimension A tot given in Table 2 by the formula: 

LAE 1 (1-x) 2 ,„ 

= r 1 - , (for x < 15 

Tivp 3 24 ' u ; v ' 

— = ^tot-\{l + X -)\ (forx>0) (16) 
Tivp 8 3 

These formulae are in excellent agreement with recent numerical renormalization group 
results of K.Ingersent and B.Jones obtained by the numerical renormalization group method 
||. They correspond to a one-parameter family of (non- Fermi liquid) fixed points (when 
symmetry between odd and even channels is preserved) which can be explored by varying 
either J or J + — J m . One of the endpoints of this line is the strong ferromagnetic fixed point 
of Ref. 0, obtained for I — > —00 (corresponding to x = 1). At this point, the impurity 
spins bind into a triplet state, and the spectrum of the 4-channel spin-1 Kondo problem is 
found, as conjectured in Ref. ||. The A(2,2) operator associated with each state at this 
fixed-point is also displayed in Table 2. The marginal line includes the point corresponding 
to two decoupled Kondo problems (x = 0) and extends over to a region of antiferromagnetic 
coupling. As discussed below, we expect it to end at x — a/6 — 3 ~ —.55 on this side. 

The above fusion principle can also be applied to find the dimensional change of any 
operator O for a non-zero x (O transforms as UgOU^), and hence the operator content at 
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a given fixed point along the line. Spin correlations can be obtained from the identification 
± S% oc de ±l ® a + h.c, which leads to a singular behaviour of the uniform susceptibil- 
ity Ximp oc T~ 9 ^ on the ferromagnetic side x > with a continuously varying exponent 
9(x) = x(2 — x)/3. Similarly, the staggered susceptibility behaves as Xst — T~ 9 (~ x ^ on the 
antiferromagnetic side x < 0. Xst (resp. Ximp) is finite for x > (resp. x < 0). The 
low-temperature behaviour of the specific heat is governed by the leading irrelevant pertur- 
bations compatible with all symmetries of the model. At the decoupled impurities point, one 
has two such operators, of dimension 3/2 which read (Si ± S^M^oo ± Je)- m the bosonised 
form, they involve de~ tyJ ^® -\-h.c and d + e~ l<i '^J2a^aJ^a + h-c, where £ a is the adjoint operator 
of the SU 4(2) algebra. For a non-zero x, they are changed into ^g-^v^+^/v^)* _|_ ^ c anc [ 
rf+e -1 *^ 1-21 )/"^ Y^ a iaJa + h.c, respectively. Hence the leading contribution to C/T is due to 
the latter for x > 0, and to the former for x < 0. This leads to C/T ~ T~ e ^ and to a 
universal (x-dependent) Wilson ratio for x > 0, while a different behaviour C/T ~ jp-^w 
with a{x) = — x(Q + x)/3 is found for x < 0. Note that this yields a = at x = x c = a/6 — 3, 
corresponding to the first operator above becoming marginal. We expect that, for x < x c , 
the system flows to the 'strong- antiferromagnetic' (Fermi liquid) fixed point where the two 
impurities bind into a singlet state, as found in ||. 

Finally, we report on some results on the effect of a non-zero coupling J_ (even-odd asym- 
metry). It can be shown that it gives rise to a relevant perturbation on the ferromagnetic 
side (x > 0), thus destabilizing the marginal line in favour of a spin-1, 2-channel Fermi liquid 
fixed point, as found in 0. For antiferromagnetic RKKY interactions {x < 0), the leading 
operators generated are marginal, and the line extends into a surface of non-Fermi liquid 
fixed points with continuously varying properties, in agreement with recent NRG findings 
These will be analysed in detail in a forthcoming work ||17|| , together with the effect of 
particle-hole symmetry breaking. 
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Table Captions 

Table 1 

Primary operators of the A(2, 2) algebra. (A) labels an operator of dimension A. The 
second line displays the bosonic representation of the operator (when it exists). NS and R 
stand for the Neveu-Schwarz and Ramond sectors of the algebra. 
Table 2 

Finite-size spectrum of low-lying states, j is the total spin quantum number, j\ (resp. 
j 2 ) is the SU (2) fi avour quantum number for I — 1 (resp. 1 = 2), and Qi (resp. Q 2 ) is the 
charge. The sixth column displays the A(2, 2) operator associated with each eigenstate at 
the decoupled impurities fixed point (/ = 0, i.e x = 0), whereas the eigth column displays 
the corresponding operator at the strong ferromagnetic fixed point (/ = — oo i.e x = 1). The 
degeneracy of each state is displayed in the last column, while A to t is the total conformal 
dimension at arbitrary x entering the formula for the excitation energy given in the text 
(eq-(E). 
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